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Abstract 

We derive two new solutions in terms of elliptic functions, one for the dark and 
one for the bright soliton regime, for the semi-discrete cubic nonlinear Schrodinger 
equation of Ablowitz and Ladik. When considered in the complex plane, these 
two solutions are identical. In the continuum limit, they reduce to known elliptic 
function solutions. In the long wave limit, the dark one reduces to the collision of 
two discrete dark solitons, and the bright one to a discrete breather. 

Keywords: discrete nonlinear Schrodinger equation, Ablowitz and Ladik model, elliptic 
function solutions. 
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1 Introduction 

Discrete versions of the nonlinear Schrodinger (NLS) systems have received tremen- 
dous attention recently, as this family of evolution equations has been demonstrated to be 
widely applicable to many physical disciplines. Although straightforward discretizations 
of NLS using second order central difference are usually nonintegrable, these models have 
served as an illustration of the relationship between disorder and nonlinearity. From the 
perspective of lattice dynamics, such studies are relevant in fields like biology, condensed 
matter physics, fiber optics, and material science [El UJ3 EE] ■ As examples for discussion, 
the Raman scattering spectra of an electronic material, the irreversible delocalizing tran- 
sition of Bose-Einstein condensates trapped in two dimensional optical lattices and optical 
fibers, instabilities in coupled arrays of waveguides are all relevant applications of the dis- 
crete NLS model [E|. In other applications, localized impurities in physical systems can 
introduce wave scattering phenomena, and excite peculiar modes at the impurity sites. 
Lattice models with repulsive nonlinear defects incorporating quintic nonlinearities have 
been studied [To] . 

A huge variety of approximation methods, e.g., finite difference schemes and variational 
techniques, has been employed for these discrete evolution equations of NLS type [El EE] ■ 
The Ablowitz-Ladik model |L is an important exception where substantial analytical 
rogress can be made, as explicit expressions for solitons and breathers can be found 



Discrete breathers (DBs), used loosely here to denote time periodic oscillations in a 
localized domain, have been studied intensively. DBs of the discrete NLS can interact 
with internal modes or standing wave phonons and exhibit a rich set of dynamics [T4] . 
DBs are also important in the consideration of self focusing and collapse phenomena. 

The focus of the present work is the special, integrable discretization of the NLS 
represented by the Ablowitz-Ladik (AL) model [2], 

.9A n A n+ i + A n _i — 2A n . A A * /A A x n 2 1 /i\ 
i-Qj- H ?2 H crA n A n (A n+1 + An-i) = 0, a =1, (1) 

with the notation A n = A(x,t),x = nh. Analytical advances for the AL and related 
models EU EDI I2U EE] will be desirable, both as useful, basic knowledge on discrete 
evolution equations, as well as providing insight for the critically important nonintegrable 
cases. As usual, the cases a > and a < in (JJJ) correspond to, respectively, the focusing 
(bright soliton) regime, and the defocusing (dark soliton) regime. The main contribution 
here is to show that several families of exact solutions for the continuous NLS have their 
counterparts in the discrete versions and to provide these analytic expressions. Solitary 
waves, i.e. reductions (x, t) — > x — ct, have been studied earlier in the literature, and 
hence attention is devoted here to other solutions. To be precise, we start with the 
doubly periodic (periodic in both x and t) solutions for the continuous NLS which can be 
expressed as rational functions of elliptic or theta functions. For simplicity, we shall call 
them bi- elliptic solutions. We choose, as an illustrative example, the solution |2~2~| \T\ lit*)] 



A 



cn(st, ki) + i\J\ + k[ dn(rx, k) sn(st, k\ 
dn(si, k%) + y/1 + k~i dn(rx, k) 



2r 2 l2 2fci 2r 2 



-iClt 



(2) 



n = z — 7-,k 



l + ki 1 + ki 1 + ki 



2 



in which the four parameters k, k\ (the distinct moduli of the Jacobi elliptic functions) 
and r, s (the wave numbers or periods in x and t) are constrained by the two indicated 
relations. In theta functions notations, an equivalent representation is 



.4 



a9^(0,r)9 A (0,r)9 2 (0,T 1 ) 

04(0,7-1) 

a 2 k 4 (0,r)+£ 4 4 (0,T)l, u 



9 2 (wt, Ti)9i(ax, r) + iOiiut, Ti)9 3 (ax, r) 
9$(ujt, Ti)9i(ax, r) + 9 A (uit, Ti)9 3 (ax, r) 



-ifit 



2a 



,9i(0,r)9l(0,T) 3 2 (O,n) 9i(0,r) + 9\(0, 



T 



(3) 



91(0, Tl ) 



9j(0, n) 2^1(0, r)fl 4 2 (0,r) 



The representation in theta functions is more symmetric and the transformation in wave 
numbers follows from the classical theories of such functions OH] • The expressions 
(J2J) and (J3J) solve the continuous NLS in the dark soliton regime, 



iA t + A xx - 2A 2 A* = 0. 



In the bright soliton regime, 



iA t + A xx + 2A 2 A* = 0, 
the corresponding solution is [221 El ED] 



(4) 



(5) 



.4 

n 



r 

V2 



(1 + ki) 1 ^ 2 dn(st, fci) cn(rx, k) + ik 1 ^ 2 sn(st, ki) 



-r 2 k u k 2 



1 — kY 2 (l + k\) x l 2 cn(st, ki) cn(rx, k) 

, 1 - kl 2 



(6) 



The method to obtain the above solutions (j2J), © is either a special assumption (221 El 
or the Hirota bilinear method JU]- in the long wave limit, the dark solution (J2J) reduces 
to the collision of two dark solitons, while the bright one © reduces to a breather. We 
shall show that both (J2J) and © have counterparts in discrete versions of NLS, or more 
precisely in the AL model (0). However, the change from continuous to discrete evolution 
equations is not completely straightforward. Even though the angular frequencies (in t) 
are still related to the wave number (in x), the constraint on the moduli of the elliptic 
functions will involve the wave number as well, see Sections 121 and El 

The paper is organized as follows. In Section [21 we show that the defocusing discrete 
NLS admits a bi-elliptic solution whose continuum limit is (J2J). In Section El the similar 
solution is presented in the focusing case. Finally, in an Appendix, we give in a quite 
symmetric notation the continuous and discrete analytic expression which unifies the 
focusing and defocusing bi-elliptic solutions in the complex plane of x and t. 



2 Discrete defocusing NLS 



A solution for © in the defocusing regime (a — —1) is 

;i - k 2 ) 1/A cn(st, ki) + i(l - kl) l/4: dn(rnh, k) sn(st, k x ) 



A, 



«ov k 



;i - k 2 ) 1 / 4 dn(st, fci) + (1 - kl) 1 / 4 dn(rnh, k) 



fci(l - k 2 ) 1 ' 2 sn 2 (rh, k) 
h 2 (l - kl) 1 ' 2 dn(rh, k) ' 
-k 2 \ 1/2 _ k 2 



Q = 2^, s = 2^ 



1-kl 



1 + dn(rh, k) ' 



-int 



(7) 



in which two parameters are arbitrary, e.g. r, k. The field A n is doubly periodic in both 
x = nh and t, see Fig. ^ A remark on the derivation is in order. It is not clear whether 
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a first order transformation described earlier in the literature [22 E] will succeed in the 
discrete version The Hirota operator relevant to discrete evolution equations will 
typically involve hyperbolic functions. Whether such an operator can be profitably applied 
here will be left for future studies. This solution (J2J) is derived here by direct application 
of identities between theta functions [SI EH]- The fact that the mere replacement of x 
by nh in (J2J) still yields a solution of the discrete NLS JTJ just reflects the remarkable 
integrability properties of ([TJ). 

The continuum limit. Using the expansions at the origin 



sn(z, k) = z + 0(z 3 ), cn(z, k) 



2^2 



Z - + 0{z% dn(z,*) = l-*y 



0{z 4 



letting nh = x and taking the limit h — > 0, the solution (J7J) of the discrete NLS goes to 
the solution (J2J) of the continuous NLS. 

The constraint on k, ki,rh in (J7J), which in the limit h — > reduces to the constraint 
on k, ki in (J2J), defines rh as an elliptic integral of (k, ki), and it admits real solutions if 
and only if k% < k. 

The long wave limit. In the limit k\ — > 1, the solution (J2J) of the continuous defocusing 
NLS (JU) reduces to the collision of two dark solitons. The quantity 
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(1 - k 2 ) 1 ' 2 

(\-k 2 yi 2 ' 



has for limit 



sech rh 
hm = : — -. 

fc->i,fci^i 1 + sechr/i 



Therefore the long wave limit of (JZj) is the collision of two dark solitons, 



A„ — A 



D~ 



S sech st + i tanh st sech rnh —i£} D t 



A 2 
D 



1 — sech rh 
h 2 



5 sech st + sech rnh 



n 



D 



2A 2 D , s 



2A 2 D , 



in which r is arbitrary. 



(9) 



(10) 



3 Discrete focusing NLS 

For with (7=1, starting from the known solution for the continuous case 

Oziujt, Ti)9 2 (ax, t) + i9i(ut, Ti)9^(ax, r) 



A = ao 



O^iujt, r^O^ax, t) — 8 2 (ujt, Ti)62(ax, r) 



(12) 



we make in the discrete case an Ansatz consistent with that solution. To satisfy (JT|). 
repeated applications of theta and elliptic functions identities now yield the solution, 



.4, 



a 



k 1/2 dn(st, ki) cn(rnh, k) + ik\ /2 {\ - k 2 ) 1/A (l - A; 2 ) 1 / 4 sn(st, h) 



-k 1 / 2 ^ 2 cn(st, h) cn{rnh, k) + (1 - fc 2 ) 1/4 (l - fc i) 1/4 
fc(l~- A; 2 ) 1 / 2 sn 2 {rh,k) 



a5 ~ /i 2 (l-A; 2 )V2 cn(r/i,A:j ' ^~ ^ 



o- 



2< 



(13) 



k 2 + fcfci 



1 - k 



2\ 1/2 



1 + cn(r/i, fc) 



in which r and k are arbitrary. We have reverted back to notations of the Jacobi elliptic 
functions in (|T3*j) as they are more compact. 

The continuum limit. Letting nh = x and taking the limit h — > 0, we recover the 
solution (jnj) through the following steps. We first expand for small h the dispersion 
relation in (fT3|) . which yields 



limk 2 = ^-^, (14) 
then, using this value of k 2 , the small h expansion of the other relations in ()13|) leads to 

lim al = — , limfi = — r 2 fci, lims = r 2 . (15) 

Eq. (|T3j) then reduces to ©• 

The long wave limit. It is of interest to take the long wave limit k\ — > 1, k — > since it 
should yield the discrete breather. We make use of the following limits of Jacobi elliptic 
functions as k\ — > 1, k — > 0, 

dn(z, k\) ~ sechz, dn(z, k) ~ 1, 

cn(z, ki) ~ sech z, cn(z, k) ~ cos z, (16) 
sn(z, ki) ~ tanh z, sn(z, k) ~ sin z. 

The two limits k\ — > 1 and — > are not independent since /i must be kept nonzero to 
get the discrete breather, and the dispersion relation in ()13|) requires 



kkAl-k 2 ) 1 ' 2 k 



;i7) 



1 + cosr/i (i-jfe^Va (1-^)1/2- 

On applying (fTTj) to (fl3~J). we obtain 

i- o i- 1-cosr/i 2 1-cosr/i 

hmiZ = - Inns = -2—- — , hma = — — . (18) 

h z cosrh h z cosrh 

The reality of ao implies 

rh 1 r- r/i 

cosrh > and cos — > — = and v 2 cos sech st cosrnh ^ 0. (19) 

2 ~ v/2 2 ^ v ; 

Applying these limits to (fE?j) . we have 

sechst cosrnh + i\/2cos(rh/2) tanhst 



^4. n — (2q _ 

sech si cos rnh + a/2 cos(r/i/2) 
2 1 -cosrh 2 2 



(20) 



an expression depending on the arbitrary parameter r, which does represent the discrete 
breather as expected. This discrete breather is localized in time, see Fig. |2l as opposed 
to another discrete breather localized in space jHJEl- 
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4 Conclusion 



A class of solutions which are doubly periodic separately in space and time has been ob- 
tained explicitly for the integrable discrete nonlinear Schrodinger equation (NLS), namely 
the Ablowitz-Ladik model. Both the focusing and defocusing regimes have been treated. 
These solutions, expressed as rational functions of the Jacobi elliptic functions, can be 
unified for both regimes into a single complex expression, using a symmetric notation 
introduced by Halphen (Appendices A and B). 

This single complex expression is a particular instance of a quasiperiodic solution 
(QPS) involving four periods in the complex plane. The mathematical problem of finding 
all the QPS of the Ablowitz-Ladik equation has been solved by various authors [7| 14^ 1231, 
and even beautifully extended to the whole Ablowitz-Ladik hierarchy by Vekslerchik [2E] . 
However, the resulting formulae are uneasy to handle for a practical use. Our solutions, 
which are derived by quite elementary algebra, are evidently included in this abstract 
mathematical framework, but they are simple closed form expressions, which involve 
two distinct elliptic functions (genus one) in the respective variables nh and t. They are 
the natural discretization of the solution of Akhmediev and Ankiewicz 

Current research efforts are made to possibly find similar solutions for nonintegrable 
discretizations of the NLS frequently encountered in physics, as well as integrable versions 
of the higher order NLS. 

An important continuous, integrable higher order NLS incorporating third order dis- 
persion and cubic nonlinearity is the Hirota equation [TBI. A finite difference version has 
been proposed and N-soliton solutions have been given [2*1] . An important link between 
the theoretical works on discrete evolution equations and nonlinear lattices exists, as in- 
trinsic localized modes in a nonlinear lattice with a hard quartic nonlinearity are governed 
by this discrete Hirota equation [TZj. An exciting new direction for future work is thus 
to examine if the doubly periodic or breather-types solutions found here can be extended 
to discrete higher order NLS. Extension and connections with nonlinear lattice dynamics 
can then be made. 
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Appendix A. Results in complex, symmetric notation 



In the complex plane of the variables x and t, the dark and bright soliton solutions, 
Eqs. (J2J) and ©, are not distinct and their common expression can be made invariant 
under any permutation of three carefully defined entire or elliptic functions, just like the 
Weierstrass function p is invariant under any permutation of the three zeros (ei, e 2 , e 3 ) of 
p' . All the solutions previously given in the text are here reexpressed in such a form. The 
relevant mathematical formulae, introduced by Halphen [T2*] . can be found in Appendix 
B. 

For convenience, let us first denote the continuous NLS, whether in the dark regime 
P|) or in the bright regime (jSj), as 



iA t + pA xx + qA A* = 0, p, q e TZ, 
and the parameters of the two elliptic functions as 

f p'(tf = 4(p(t) - E 1 )(p(t) - E 2 )(p(t) - E 3 ) = Ap\t) - G 2 p(t) - G 3 , 
1 p'O) 2 = 4(p(ar) - e 1 )(p(x) - e 2 )(p(x) - e 3 ) = Ap 3 (x) - g 2 p(x) - g 3 . 



(21) 
(22) 



Finally, let (a,/3, 7) and (a, b, c) be two independent permutations of (1,2,3). 

The bi-elliptic dark or bright solution, Eqs. or is then represented by the 
symmetric expression 



, _ aih a (t) + ia 2 hb(x) h 7 (£) _ int _ aicr(x)a a (t) + ia 2 a b {x)a^{t) _ int 

A — 7 — n 6 — 7 — \ 7 — x T\ 6 1 



hp{t) + «4 hfe(x) 
tt 2 _ 2p 

q q 



a 



[x)ap(t) + a i a b {x)a{t) 



-qaia 2 , a 2 = 2pVt, 



(23) 



Ep — E~ = f2 , 3e b 



f2 En, — E~ 



P Es-E, 



3e b 



in which the two arbitrary constants are, for instance, (e c , Q) . The expression (J23j) is 
invariant under any permutation of the indice s (a , /3, 7) and (a,b,c) independently. 
In the discrete case, the solutions (J7|) and (|T3j) can both be represented as, 



A, 



ai h a (t) + ia 2 h. b {nh) h 7 (£) aia{nh)a a {t) + ia 2 a b (nh)a^(t) 



-int 







a\ = 




q 


%~ 


Ej 


E a - 


E^ 


Ep — 


E^ 



hp(t) + a 4 hb(nh) 
2 2p 



Vl\ 2 



qh 2 h a (h) hc(h) 
K(h)h c (h)-ht(h) n 



a{nh)ap{t) + a4a b (nh)a(t) 

2 2pQ 
-qaia 2 , a 4 



(24) 



h 2 h a (h) hc{h) 



9e 2 



P 



2 ■ 



4(h a (h)hc(h)-h 2 b (h)) 

in which the two arbitrary constants are, for instance, (e c ,f2). In the continuum limit 
h — > 0, according to (|3*K|) . one has 



lim h 2 ha(h) hc(h) = 1, ha(h) h c (h) - h 2 b (h) 

h— >0 



3e b 
2 



e c ) 2 y + 0(/* 4 ), (25) 



so ()23j) goes straightforwardly to f|23j) . 

One might wonder whether the trigonometric degeneracies (vanishing of the elliptic 
discriminant) occur independently in x and t or simultaneously. In the discrete case, the 
discriminants are 



g\ - 21 g\ = 16 [(e b - e c )(e c - e a )(e a - e b )\ 



Gl 



21G\ 



256(e a - e b ) 2 (e b - e c 



9e 2 + 4(h a (h) hc(h)-hl(h) 



Q"2 ^4 

af 



When the last factor in G% — 27Go vanishes, i.e. 



9e 2 b + 4 (ha(h) hc(h) - h 2 b (h)) = 



(26) 



by elimination of h a (^)h c (^) and hj,(/i) with the rules of derivation (J32"j) and ({37}), this 
implies (e a — e c ) h&(/i) = 0. Therefore, the trigonometric degeneracy occurs simultaneously 
for x and t. By the continuum limit, this is also true in the continuous case, where the 
two discriminants evaluate to 



gl 



Gl 



27G 



2 _ i ( a 2 a!\ 4 / 
\a x al) \ 



C3 



12a? 



e 3 + 



12al 



(#2 " 27<7. 



(27) 
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Appendix B. The symmetric notation of Halphen 

The notation of Jacobi (four entire functions $1,2,3,4 and twelve elliptic functions pq, with 
p and q chosen among s,c,d,n) makes the practical computations quite technical, because 
of 

1. the lack of symmetry under a permutation of the indices 1, 2, 3 of dj or of the letters 
c,d,n of the twelve functions pq, 

2. the lack of homogeneity, which implies that the arguments u, v, w involved in the 
relations between p(u), ^(f ), pq(u>) are all different, 

3. the high number (twelve) of elliptic functions thus defined. 

All these inconveniences have been removed by Halphen , but his elegant notation 
is totally absent from the handbook [Bj (the notation of Neville is different), only presented 
in a less symmetric way in p2], and scattered among different locations in [IBj, therefore 
we summarize it in the present Appendix. 

In addition to the odd entire function a of Weierstrass, Halphen defines three even 
entire functions a a ,a = 1,2,3, 12., Chap. VIII p. 253] [H §10.5 p. 391] 18, Eq. (6.2.18)] 

a a (u) = + u ) e -VaU } Va = c(Wo) Q = 1 2 3 (28) 

<J{0J a ) 

with the usual notation (HI §18.3.4] 

uji=u, u 2 = u + u', u 3 = u', 771 = 77, r]2 = V + V', Vs = V, ( 29 ) 

in which 2u, 2u/ are the two periods. Their symmetry with respect to permutations of the 
indices results from the relation [H Chap VI p. 191] [19 , §10.5 p. 391] [HI Eq. (6.9.1)] 

J p{u) -e a = ^7"T' « = 1, 2, 3, (30) 

To avoid any ambiguity with the square roots, let us define the three odd functions 

ha(u) = ——. (31) 
a(u) 

By construction, all formulae involving h a are invariant under any permutation (a, (3, 7) 
of (1,2,3). Then, the only relations needed to establish all the results of Appendix A 
are: 

the derivation formula 

l4(tt) = -h / j(«)h 7 («) > (32) 

the algebraic dependence relations 

hp(u) + ep = hj(u) + e 7 , (33) 

the addition formula 

ha ( u + v \ = h'( M )h'(^) ~ (e a -e /3 )(e a -e 7 ) 

ha(u) hp(v) h 7 (f ) + h a {v) hp(u) h 7 (u) ' 
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the Laurent expansion at the origin [T2J Chap VII p. 237] [121 Chap IX p. 304] 
, / \ 1 u . / , 2 \ « 3 %3 - 7e a g 2 5 
28^ - 225e Q ^ 3 - 105e^ 2 7 



+ 



V + 0(« J 



2 9 .3.5 2 .7 

and the degeneracy to trigonometric functions O Chap VIII (75) p. 289] 



(35) 



7T- 4 / 7T V 



^ 12' ^ 2 3 V2cj7 



#3 



27 V2w 



7T 



1 / 7TU 



a^wj = a 7 ^j = exp i - 



e/3 = e n 



7VU 



3#3 

2<? 2 : 



2cj 



7TU 



(36) 



cx Q (wj = cr^ttj cos — , a\u) = — a p [u) sin — . 

LUJ 7T LUJ 



From (J32*Jl and (jSSJ), one deduces the differential equation 

{h' a {u)) 2 = (h 2 a (u) + e a - e p )(Tol(u) + e a - e 7 ), 



(37) 



and from (J34)) the expression of the three elliptic functions h a as rational functions of p 
and p' (a characteristic property of the elliptic functions), 



(p(u, 92, 93) - e a ) 2 - (eg - eg) (e a - e 7 ) 

p'(u,g 2 ,g3) 



(3* 



The link with the asymmetric notation of Jacobi involves four relations between the 
two sets of four entire functions, [H Chap VIII, (49) p. 260] [H §10.5] 



2 (v) 

3 («) 



^(e 2 -e 3 ) 1 / 4 e-^ 2 /(2-) (Tl(M)) 
2W (e 1 -e 3 ) 1 / 4 e-^ 2 /(2-) (T2(u)? 
( ei -e 2 ) 1 / 4 e-^ 2 /(2^) (T3(u)? 




(39) 



7T 

2uj 




A = 16(ei - e 2 ) 2 (e 2 - e 3 )> 3 - e^, k 2 = m 



e 2 - e 3 
ei - e 3 : 



plus four relations (among them a scaling one) between the trio h a and the copolar trio 
ps(u), p=c,d,n [EI Chap II, (16) p. 46], 



cs(z) ds(z) ns(z) u 1 
hi(w) h 2 («) h3(«) z s/ex - e 3 ' 



(40) 
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5 Figures captions 



Figure 1: Intensity \A n \ 2 vs. x and t for the discrete defocusing bi-elliptic solution (0), 
with parameter values k = 0.5, k\ = 0.6, r = 1, h = 0.2. 



Figure 2: Intensity \A n \ 2 vs. x and t for the long wave limit (J20|) of the discrete focusing 
bi-elliptic solution, with parameter values r = 1, h = 0.2. 
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